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ABSTRACT 

We present new Spitzer Space Telescope observations of two fields in the Orion 
OBI association. We report here IRAC/MIPS observations for 115 confirmed 
members and 41 photometric candidates of the ~10 Myr 25 Orionis aggregate 
in the OB la subassociation, and 106 confirmed members and 65 photometric 
candidates of the 5 Myr region located in the OB lb subassociation. The 25 
Orionis aggregate shows a disk frequency of 6 % while the field in the OBlb 
subassociation shows a disk frequency of 13 %. Combining IRAC, MIPS and 
2MASS photometry we place stars bearing disks in several classes: stars with 
optically thick disks (class II systems), stars with an inner transitional disks 
(transitional disk candidates) and stars with "evolved disks"; the last exhibit 
smaller IRAC/MIPS excesses than class II systems. In all, we identify 1 transi- 
tional disk candidate in the 25 Orionis aggregate and 3 in the OBlb field; this 
represents ~10% of the disk bearing stars, indicating that the transitional disk 
phase can be relatively fast. We find that the frequency of disks is a function 
of the stellar mass, suggesting a maximum around stars with spectral type MO. 
Comparing the infrared excess in the IRAC bands among several stellar groups 
we find that inner disk emission decays with stellar age, showing a correlation 
with the respective disk frequencies. The disk emission at the IRAC and MIPS 
bands in several stellar groups indicates that disk dissipation takes place faster 
in the inner region of the disks. Comparison with models of irradiated accretion 
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disks, computed with several degrees of settling, suggests that the decrease in the 
overall accretion rate observed in young stellar groups is not sufficient to explain 
the weak disk emission observed in the IRAC bands for disk bearing stars with 
ages 5 Myr or older; larger degrees of dust settling are necessary to explain these 
objects. 

Subject headings: infrared: stars: formation — stars: pre-main sequence — open 
clusters and associations: individual (Orion OBI association) — protoplanetary 
systems: protoplanetary disk 



Introduction 



Observational and theoretical studies indicate that important processes in the evolution 
of protoplanetary disks take place at ages between 1 and 10 Myr. Ab out 90% of the low 
mass stars (~K5 or later) have lost their primordial disks at 5-7 Myr (e.g-. lHaisch et al.ll2001 



Hartmannl 120051 : IHernandez et al.l 120071 ). Grains grow to sizes of ~1000 km stirring up the 
leftover small object s in the disks and o r iginating the first gene ration of reprocessed dust by 



collisional cascades (IKenvon et al. 



to form in this period (jPollack et al 



2005 



Hernandez et a. 



199 



6; 



Alibert et al. 



2006). Giant planets are expected 



2004J ) . However, additional studies 



of disk population in this crucial age range are necessary to improve our knowledge and 
clarify many details about the evolution from primordial disks to planetary systems. 

OB associations are excellent laboratories for comparative studies of protoplanetary 
disk evolution, because they harbor young stellar populations (1-10 Myr) originating from 
the same giant molecular clouds, spanning a wide range of ste l lar masses, and in a vari- 



ety of evolutionary stages and environments (jBrown et al.l Il999i : ISicilia-Aguilar et al.l 12006 



Briceno et al.ll2007al : iPreibisch fc Zinneckerl 120061 ) . In particular, the Orion OBI association 
(Ori OBI), as other OB associations, s hows a well defined age sequence suggesting a large- 
scale triggered star formation scenario (jBriceno et al.ll2005l l2007al : iLee fc Chenl 120071 ). Ori 
OBI contains very youn g subgroups (ages < 1 Myr) still embedded in their natal gas (e.g., 



Orion A and B clouds; iMegeath et al.l 120061 ) . sub groups in the process of dispersing their 



natal gas (e.g., the a Orionis cluster, age ~ 3 Myr; IHernandez et al.ll2007l ) and more evolved 
populations, which have long sinc e dissipated their pro genitor molecular clouds (e.g., the 25 
Orionis aggregate, age~ 10 Myr; iBriceno et al.ll2007bl ). 



We are carrying out an optical photometric and spectroscopic survey of ~ 128 deg 2 in 
Ori OBI in order to identify the low and intermediate mass stell ar populations, and study 



the properties linked to the first stages of star and disk evolution (jBriceno et al.ll200ll . 12005 
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2007bl Jd: ICalvet et al.l l2005al : IHernandez et al.l I2005L 120061 . l2007h . In this work, we expand 



the results from the optical survey with the capabilities of the Spitzer Space Telescope at 



near and mid infrared wave 



engths to identify and characterize pro t oplanetary disks around 



young stellar objects (e.g. 


Allen et al 


2004 


; 


Meeeath et al. 


2004 




Gutermuth et al. 


2004 




Muzerolle et al. 


2004 


Hartmann et al.l 


2005|; 


Sicilia-Aeuilar et al. 


2006|; 


Hernandez et al.ll2007) 



In particular we study the near and mid infrared properties of stars in two IRAC/MIPS 
Spitzer fields encompassi ng an area of ~2.5 de g 2 . One field is located in the 7-10 Myr 25 
Orionis stellar aggregate (jBriceno et al.ll2007bl ). the most populous ~ 10 Myr stellar group 
known within 500 pc; the other is located in the Ori OB lb sub-association, in whi ch we 
have estimated an age of ~ 5 Myr (jBricefio et al.ll2005l. l2007bl; IHernandez et al.ll2006l ). Ad- 
ditional results from the 3 Myr a Orionis cluster (IHernandez et al.l 120071 ). also located in 
the Ori OBlb subassociation, allow us to cover most of the potentially crucial age range 
in protoplanetary disk evolution. In this cluster, we found 336 photometric members using 
optical and near infrared color-magnitude diagrams, about a third of this sample exhibits 
excess in the IRAC and/or MIPS bands indicating that they have disks. 

This paper is organized as follows. In §2] we present the observational data and a brief 
description about membership. We analyze the observations and describe the disk emission 
detection in 331 The main results are summarized in 



2. Observations 



2.1. Infrared photometry 

We have obtained near-infrared (NIR) and mid-infrared photometry of two regions in 
the Orion OBI associat ion using the four channels (3.6, 4.5, 5.8 & 8.0 /zm) of the InfraRed 
Array Camera (IRAC . iFazio et al.l 120041) . and the 24 \xm band of the Multiband Imaging 



Spectrometer (MIPS; iRieke et al.l l2004j ). on board the Spitzer Space Telescope. The field 
located in the 25 Orionis aggregate (hereafter "25 Orionis") covers an area of ~1.1 deg 2 
centered at RA~ 5.42 hours and DEC~1.64 deg; the other field (hereafter "OBlb") covers 
an area of ~1.4 deg 2 on the Orion OBlb sub -association centered at RA~5.52 hours and 



DEC~-1.71 deg. Dust infrared emission maps (jSchlegel et al.lll998l ) reveal that at least 90% 
of the regions covered by IRAC images in 25 Orion is and OBlb have vis ual extinctions 
smaller than Ay~0.12 and Av~0.6, respectively (see IHernandez et al.l 120061 ) . These values 
are mostly in agreem ent with the mean visual extinction calculated from individual stars in 
Briceno et all J2005I ). 



The IRAC observations were done using a standard raster map with 290" offsets, to pro- 
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vide maximum areal coverage with just a slight overlap between frames, to aid in mosaicking 
the data. Each position is composed of 3 dithers, with a single-fra me integration of 12 



seconds. The IRAC observations were processed using the IRACpro c (ISchuster et a 



package to create the final mosaics with a scale of 0.86 "/pixel (see IHernandez et al 



2006 



1) 

2006( 1. 



Point source detections were carried out individually on each IRAC channel using Phot Vis 
tool (an IDL GUI-based photometry visualization tool developed by R. Gutermuth). More 
than 20,000 sources in each field were detected in at least one Spitzer band. We extracted the 
photometry of these objects using the apphot package in IRAF, with an aperture radius of 
3". 7 and a background annulus from 3.7 to 8". 6. We adopted zero-point magnitudes for the 
standard aperture radius (12") and background annulus (12-22". 4) of 19.665, 18.928, 16.847 
and 17.391 in the [3.6], [4.5], [5.8] and [8.0] channels, respect ively. Aperture cor rections were 
made using the values described in IRAC Data Handbook (jReach et al. II2006I ). 



MIPS observations were obtained using the medium scan mode with full-array cross- 
scan overlap, resulting in a total effective exposure time per pointing of 40 seconds. The 
images were processed using the MIPS instrument team Data Analysis Tool (DAT), which 
calibrates the data and applies a distortion correct ion to each individual exposure before 
combining it into a final mosaic (IGordon et al.ll2005l ). We obtained point source photometry 
at 24 fj,m with IRAF/it daophot point spread function fitting, using an aperture size of 
about 5.7" and an aperture correction factor of 1.73 derived from the STinyTim PSF model. 
The absolute flux calibration uncertainty is less than 5%. Our final flux measurements are 
complete down to about 1 mJy in both maps (the limit flux is about 0.5 mJy). 

Figures [T] and [2] show color images combining three channels of IRAC ([3. 6], [4. 5] and 
[8.0]) for 25 Qr ionis and for Ori OBlb, respe ctively. We display the low mass spectroscopic 
members from iBriceno et al.l (120051 . l2007bl Jd) and the low mass photometric candidates se- 
lected in §2.31 the stars bearing disks studied in §3.1( and the intermedia te mass members 



i nclud ing the debris disk candidates and the Herbig Ae stars studied in IHernandez et al. 
j2006h . 



2.2. Optical photometry 



Optical (V and I Cousin) magnitudes were obtaine d from the CIDA V ariability Survey 
which is being carried out using the QUEST I camera (IBaltay et al.ll2002l ) installed on the 
Jurgen Stock Telescope (a celar aperture 1-m Schmidt telescope) at the Venezuela National 
Astronomical Observatory. The camera, an array of 4x4 CCDs, is designed to work in 
driftscan mode, which is a very efficient way to survey large areas of the sky. Each scan 
was reduced and calibrated with the standard QUEST software and the method described 
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in 



Vivas et al.l (120041 ). in which variable stars can be identified (see IBriceno et al.ll2005l ). 



2.3. Low mass members and photometric candidates 



We follow the procedures described in IHernandez et al.l (120071 ) to reject non stellar 
objects and contaminating sources using IRAC color-color and IRAC color-magnitude dia- 
grams. In brief, we select stars with [3.6] < 14.5, below this limit, the co ntamination from 
extragalactic sources is expected to be more than 50% (IFazio et al.ll2004l ). The [4.5] — [5.8] 
versus [5.8] — [8.0] and [3.6] — [5.8] versus [4.5] — [8.0] color-color diagrams were used to elimi- 
nate most of the galaxies with polycyclic aromatic hydrocarbon (PAH) emission and objects 
with strong 8 yumcontamination (Gutermuth et al. 2007, in preparation). 



Optical and 2MASS counterpart (jCutri et al.l 120031 ) for the IRAC sources were found 
using a 2" matching radius. A preliminary list of 623 objects in 25 Orionis and 918 objects in 
OBlb were created using optical-2MASS color magnitude diagrams (V versus V-Ic, V versus 
V-J and J versu s J-K) to select those objects above the zero age main sequence (ZAMS; 
Siess et al.l lioool l at the distance of each stellar group (330 pc and 440 p c for 25 Orionis 
and OBlb, respectively; IBriceno et al.ll2005l l2007bl ; IHernandez et al.ll2005l ). We rejected by 
visual inspection non-members sources, like diffuse objects, and objects with an apparent 
problem in the photometry, like close binaries, faint companion binaries, and stars on the 
image border. 

Low mass members were confirmed by Bricefio et al. ( 2005 . 2007b q) using optical spec- 
tra obtained at the Fred Lawr ence Whipple Observa tory with the 1.5 m telescope equipped 
with the FAST spectrograph (IFabricant et al.lll998l ). and with the 6.5 m MMT telescope 



equip ped with the Hectospec and Hectochelle multifiber spectrographs (IFabricant et al. 
20051 ). Additional low resolution spectra were obtained at the Kitt Peak National Ob- 



servatory using the Hydra multiobject spectrograph on the 3.5 m WIYN telescope. Low 
mass members of 25 Orionis and Ori OBlb can be identified by the p resence of Li I A670 7 
in absorption and the equivalent width of the Ha line in emission (e.g. IBriceno et al.l 119981 ). 



The strength of the Ha line is used to separate accreting stars, the Classical T Tauri stars 
(CTTS), from non-accreting stars, the Weak T Tauri stars (WTTS). To distinguish b etwee n 
the two type of objects we followed the prescription adopted by IWhite fc Basril (120031 ). 
Spectral types were derived using the SPTCLASS code@, which includes several molecular 
features (like TiO and VO bands), characteristic of low mass stars. Some stars in our sam- 
ple have additional spectroscopic confirmation using radial velocity distributions. We refer 



http: / /www. astro. lsa.umich.edu/hernandj/SPTclass/sptclass/html 
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the reader to IBriceno et al.l (120051 . l2007bl id) for more details about spectroscopic member- 
ship confirmation of stars belonging to the Ori OBI association. Tables [1] and [2] show the 
IRAC and MIPS photometry for 115 confirmed members of 25 Orionis and 106 confirmed 
members of Ori OBlb. Column 1 shows the intern al running identifi cation number in each 
sample, columns 2 shows the 2MASS designation (ICutri et al.l 120031 ); columns 3 and 4 are 
the stellar coordinates; columns 5, 6, 7 and 8 give IRAC magnitudes in the bands [3.6], [4.5], 
[5.8] and [8.0], respectively; column 9 gives MIPS (24 fan) magnitudes; column 10 gives the 
reference for the optical spectroscopic and photometric data ; the last column gives the disk 
classification based on the IRAC and MIPS analysis ( §3.1)) . 

Figure [3] shows optical-2MASS color magnitude diagrams (CMDs), V vs V-I and V ver- 
sus V-J, illustrating the selection of the photometric candidates o f 25 Orionis (upp er panel) 
and Ori OBlb ( lower panels). Solid lines represent the ZAMS (jSiess et al.l 120001) for eac h 



stellar group. Dashed lines represent the 10 Myr and 5 Myr isochrones (jSiess et al.l 120001 ) . 



Since, it is well 
ate mass stars ( 



mown that theo retical, non-birthline isochrones do not so well for intermedi- 
Hartmarml 2003 ) and the opacity tables at low temperatures are incomplete 



( iLyra et al.l 120061 ; iBaraffe et al.lll998l ). these theoretical isochrones are plotted as reference 



and do not affect the selection of the ph otometric candidates. We a s sumed a distance of 33 
pc for 25 Orionis and 440 pc for OBlb (IBriceno et al.ll2005l . l2007bl ; Hernandez et al.l 120051 ) . 
The distribution of confirmed members (open circles) can be roughly traced as a straight line 
in the CMD diagrams, with a larger spread in Ori OBlb due to the larger visual extinction 
in this region in comparison with 25 Orionis. We calculated the median colors (V-I and V-J) 
for confirmed members (open circles) within 1.0 mag bins in the visual band. We used the 
differences between the observed colors and the expected colors (the median; represented 
by long dashed lines) to calculate the standard deviation (a) of the member samples. We 
selected as photometric candidate s (open squares) stars located in the region described by 



the 2.5ct limits in the CMDs (see Hernandez et al.ll2007f ). Tables |3J and H show the IRAC, 



MIPS and optical photometry for the 41 photometric candidates found in 25 Orionis and 
for the 65 candidates found in OBlb, respectively. The information shown in these tables in 
columns from 1 to 9 is the same as in Table [H Columns 10 and 11 show the V magnitude 
and the color V-I from the CIDA survey; the last column gives the disk classification based 
on the IRAC and MIPS analysis (333]). 
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Results 



3.1. Disk diagnostics 



Figures H] and show three diagrams used to identify and roughly characterize the stars 
bearing disks in 25 Orionis and OBlb, respectively. The top panels show the SED slope, 
determined from the [3.6] — [8.0] color (IRAC SED slope), versus the [8.0] magnitude for 
members (open circles) and photometric candidates (open squares). The photospheric level is 
described by the upper soli d line, which is calculat ed using the photometric errors propagated 
from the [3.6] — [8.0] color ([Hernandez et al.l 120071 ). Stars with excess emission at 8 /zm can 
be identified in this d iagram. For comparison, IR AC SED slope histogra ms for Taurus 
( jHartmann et al.ll2005l ) and for the a Orionis cluster ([Hernandez et al.l 120071 ) are displayed. 
In general, disk bearing stars in Taurus (solid histogram) are located in a well defined regio n 
(which we call the class II region) with an IRAC SED slope > -1.8 (see lLada et al.ll2006l ): 
this limit ( dashed lines) is used to identify objects with optically thick disks in which the 
inner disk emission has not been affected significantly by evolutionary processes. In contrast, 
15% of the disk bearing stars in the a Orionis cluster exhibit smaller IRAC excesses (dashed 
histogram) sugg esting a reduction in di sk photosphere height, possibly due to dust growth 
and/or settling (IHernandez et al.l 120071 ) . The bottom left panel shows the K— [24] versus 
V-J color-color diagram, in which we identify members (big open circles) and photometric 
candidates (big open squares) with 24 /xm infrared emission above th e photospheric level 



(solid lines) indicating that disks are present around these objects (e.g.; iGorlova et al.ll2006 
Hernandez et al.ll2006l . l2007l ). In this panel, we display the K— [24] color distribution for stars 
bearing disks in the a Orionis cluster with an IRAC SED slope > -1.8 (which represents a disk 
population similar to those found in Taurus) and we use this histogram to identify stars with 
K— [24] color characteristic of stars with optically thick disks ( K— [24]>3.5, class II region). 
In the upper panel and in the bottom left panel, we define the "evolved disk region" between 
the class II region and the photospheric region. The bottom right panel shows the IRAC 
color-c olor diagram, in which we identify sta r s with excess emission in the IRAC bands 



ag_. 



Allen et al. 



2004 



2006; Hernandez et al. 



Megeath et al.l 120041 ; lHartmann et al.l 120051 ; ISicilia-Aguilar et al. 



20071 ). The dashed box displays the color s predicted for CTTS of 



different accretion rates by the models of iD'Alessio et al.l (j2005bl ). In ge neral, the IRAC 



color s observed for disk bearing stars in Taurus are located in this region (jHartmann et al. 
20051 : ISicilia-Aguilar et aDbood ). 



In the top panel of Figure HI we identify six members and two photometric candidates 
located on the IRAC class II region; most of them are located near the class II limit possibly 
indicating that these objects have begun the process of clearing the inner primordial disk. 
Two members have very small IRAC excesses just above the photospheric region. These 
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objects, also located between the photospheric and the CTTS regions in the IRAC color-color 
diagram, have no MIPS detections and therefore it is not clear if the small excess observed 
at 8 fim originates from PAH background contamination, by an unresolved companion, or by 
disks present around these stars (flagged as "disk [8]?" in Tabled]). Of particular interest are 
the member la_1121 and the photometric candidate la_1626 which are located between the 
photospheric and the class II region in the V-J versus K— [24] diagram indicating that the 
outer disks around these objects are in a more evolved stage. Moreover these stars are also 
located on the photospheric region in the IRAC color-color diagram and in the IRAC SED 
slope diagram indicating that the inner disk has already dissipated and no disk emission 
can be detected at wavelength <8 ^m. The star la_1626 also has a very small excess at 
24 /im ( ~ 2 a above the photospheric level) indicating that the presence of a disk around 
this object is not yet conclusive. Overall, in 25 Orionis we identify 7 stars with disks in the 
member sample (disk frequency 6.1±2.3%), and 3 in the photometric candidate sample (disk 
frequency 7.3±4.2%). 

Similarly, in the top panel of Figure [5] we identify 13 members and 4 photometric can- 
didates in Ori OBlb that show IRAC and MIPS excesses; five of these objects are located 
between the class II region and the photospheric region. Eight members and one photo- 
metric candidate with no MIPS detection are located in the evolved disk region (flagged as 
"disk[8]?" in Table [2] and 0J . The existence of disks around these objects needs additional 
confirmation since they could be below the MIPS detection limit, or could be contaminated 
by PAH background emission (in Figure [2] it can be clearly seen that the sky background 
emission at 8 /im is very patchy, and significant at some locations). In general, the range 
of infrared excesses at 24 \xm in OBlb is similar to that of the optically thick disks in the 
a Orionis cluster. Only one disk bearing star (6%), the star lb_337, has 24 /im excess be- 
low the class II limit while 6 stars (39%) have 8 /im excess below this limit. This suggests 
a more rapid decrea s e in dust emission in the inner disk, in agreement with results from 



Sicilia-Aguilar et al.l (120061 ) in the Cepheus OB2 association. The member lb_337, located 



in the evolved disk region in the V-J versus K— [24] diagram, does not have excess in the 
IRAC bands. Overall, in Ori OBlb we identify 14 stars with disks in the member sample 
(disk frequency 13.1±3.5%) and the 4 disk systems in the photometric candidate sample 
(disk frequency 6.2±3.1%). 

Figure [6] displays the distribution of the disk bearing stars in a SED slope space diagram 
for 25 Orionis (left panel) and OBlb (right panel). The vertical axis is the SED slope 
calculated from the K— [5.8] color and the horizontal axis is the SED slope calculated from 
the K— [24] color. The dashed areas define the photospheric level calculated with the STAR- 
PET Spitzer tool for stars K5 or later. By compa rison the low mass stars bearing disks of the 
3 Myr a Orionis cluster ([Hernandez et al.l 120071 ) are also plotted. Error bars represent the 
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quarti les of disk bearing stars in Taurus, calculated from the K— [5.8] color in lHartmann et al. 



(2005) and from the K— [24] color estimated from the me dian SED slope in iFurlan et al. 



( 120061 ). in the a Orionis cluster from IHernandez et al.l (120071 ). in 25 Orionis and in Ori OBlb 
from this work. An overall decrease in the infrared excess is observed from the 1 Myr old 
stars in Taurus to the 5 and 10 Myr old stars studied in this work; the disk population of 
the 3 Myr a Orionis cluster represents an intermediate stage in this evolution. Moreover, it 
is apparent that the decrease is larger at 5.8 /im than at 24 /im, indicating that evolution 
processes occur faster in the inner region of the disk. 

In order to characterize the stars bearing disks in 25 Orionis and OBlb, we identify 
several regions in Figure [6] defined by the dotted lines ("class II region", "evolved disk 
region" and "transitional disk region"). The horizontal dotted line represents the lower 
quartile of the a Orionis cluster. Above this line ~96% of the stars bearing optically thick 
disks in Taurus are also located, indicating a limit where the inner disk emission has not been 
affected significantly by evolutionary process. We can identify the class II objects as stars 
located above this line. In general, the class II objects identified using Figure [6] are located in 
the class II region in Figures HI and [51 Disk bearing stars below the dotted lines have decreased 
the disk infrared emission at 5.8 //m due to a decrease in the irradiation surface of the inner 
disks, and so they are in an stage where processes for inner disk dissipation have begun. The 
vertical dotted line represents the lower quartile of the stars bearing disks in the a Orionis 
cluster (the lower quartile of Taurus is rightward from this line). Using this limit, the stars 
located below the dotted lines could be sub-grouped based on their disk emission at 24 fxm: 
"evolved disks objects" ( SED slope K— [24] < -1.2), in which we see an overall decrease in 
the disk emissi on in the IRAC and MIPS bands, indic ating; similar evolution in the inner 
and outer disk (ILada et al.ll2006l ; IHernandez et al.ll2007l ); and "transitional disk candidates" 
( SED slope K— [24] > -1.2), whi ch have an inner opt ically thin disk region, combined with an 
outer, optically thick dis k (e.g.; ICalvet et al.ll2005bh . As re ference, we plotted 3 transitional 
disk stars, Coku Tau/4 



D'A 



essio et al 



2005a|), TW Hya (ICalvet et al.l 120021 : lUchida et al. 



2005bh . which occupy the region defined for the "transitional 



2004|), GM Aur (ICalvet et al. 
disk candidates". In brief, we identify 5 stars with optically thick disks (class II objects), 
one transitional disk candidate and 4 evolved disk objects in 25 Orionis. We also identify 

10 class II objects, 3 transitional disk candidates and 5 evolved disk objects in OBlb. In 
spite of the evolved disks and transitional disks objects being a subsequent stage from class 

11 objects, it is not clear if transitional disk objects are a pre-stage of evolved disks or each 
stage represents an independent stage from class II objects. 

Figure [7] shows SEDs for selected stars in our samples illustrating the disk classification 
based on Figure El The first row of SEDs shows stars with optically thick disks (CII) located 
above the dotted line in Figure [61 The second row shows transitional disk candidates (TD) 
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located right and below the dotted lines in Figure [6j Finally, last two rows of panels show 
stars with evolved disks located left and below the dotted lines in Figure [6j 



3.2. Models 



We have calculated SEP slopes for models of irradiated accretion disks including dust 
settling from ID'Alessio et al.l (120061 ). In these models the disk is assumed to be steadily 
accreting at rates of M = le-9, le-8, and le-7 M©/yr, onto a star with m ass of 0.6 an d 
luminosity of 1.2 L , which corresponds to a K7 star with age of 1 Myr (ISiess et al.l I200CH ) . 
Dust settling was included using two populations of grains (big and small grains) having 
different spatial distributions, with the larger grains concentrated toward the midplane. The 
small grains located in the upper layers have different depletions given by the e parameter 
(with values= 1, 0.1, 0.01, 0.001), which is the ratio of the dust to gas mass ratio of sma ll 



grains relative to the the standard dust to gas mass ratio {(small / Cstd] ID'Alessio et al.ll2006l ). 
The inner wall of the disk, located at the dust destruction radius, was settled self-consistently 
with the same degree of depletion used in the outer disk. 

Figure [8] shows the theoretical SED slopes derived from the colors, [3.6] — [8.0], K— [5.8] 
and K— [24] versus the degree of settling represented by e. SED slopes were calculated 
convoluting the theoretical SED with the transmission curves of the respective filters. We 
plot two inclination angles along the line of sight , 30 deg (left panels) and 60 deg (middle 
panels); this range in angles represents 40% of probability of observation. Accretion rates are 
indicated for the different curves plotted in each panel. The slope has a strong dependency 
on M, showing flatter slopes for the fastest accretors; the smallest variation in disk emission 
with M is observed for the slope K— [24] of disks without settling (e=l). In general, models 
with M = 10~ 9 M Q /yr show a stronger dependence with dust settling than models for large 
accretion rates. 

By comparison, we plotted in the right panels of Figure [H] the quartiles observed for 
disk bearing stars in Taurus, in the o Orionis cluster, in Ori OB lb and in 25 Orionis. 
The range of disk emission observed in Taurus ( 1-2 Myr) can be explai ned by the models , 
indicating optically thick disks systems w i th several degrees o f settling ( Furlan et al.lbooij ) 
and accretion rates (IHartmann et al.lll998l ; ICalvet et al.ll2005al ). Most of the stars with disks 
in the a Orionis cluster (~80%) can be explained by the theoretical SED slopes but with 
small accretion rates or/and higher degree of dust settling than in Taurus. Approximately 
half of the disks observed in 25 Orionis (Figure H|) and OBlb (Figure |5|) require models 
with lower accretion rates (M < 10~ 9 M Q /yr) or/and large degree of settling (e<0.001) to 
explain the weak disk emissions observed at [3.6] — [8.0] and K— [5.8]. However, ~75% of disk 
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bearing stars in 25 Orionis and OBlb have disk emissions at K— [24] in agreement with the 
SED slopes predicted by the models, supporting the scenario where the inner disk evolves 
faster than the outer disk. 



3.3. Disk frequencies 



In Figures HI [5] and [6] we identify 7 members bearing disks in 25 Orionis and 14 in 
Ori OBlb, indicating disks frequencies in the member samples of 6.1±2.3% and 13.1±3.5%, 
respectively. These frequencies include objects with 24 /im excess as disk bearing stars. In 
§3.1} we also identify 2 members of 25 Orionis and 8 members of Ori OBlb that exhibit 
IRAC excesses but have no MIPS detections; if we add these stars as members with disks, 
the disk frequencies increase to 7.8±2.6% in 25 Orionis and to 20.6±4.4% in Ori OBlb. 
These later values are in better agreement with the disk frequencies calculated for the low 
mass stars in the Ori OBla (7 ± 3 %) and Ori OBlb (17 ± 4 %) sub associations, using 
the excess emission (2a above the expected photospheric level) from the 2MASS color H-K 
( [Hernandez et al.ll2005l ). The lower disk frequencies derived using the 24 ^xm excess could 
indicate a possible observational bias pro duced by the flux limit of MIPS observations (~ 0.5 
mJy). Assuming the distances and ages (jSiess et al.ll2000l ) plotted in Figure [31 with a visual 
extinction of Av=0.12 mag for 25 Orionis and Av=0.6 mag for OBlb, we cannot expect to 
detect disks around a O.6M star in the 25 Orionis group if Ep4] <2.5, and if E[ 2 4] <4.5 for 
objects in Ori OBlb; where the excess ratio, E[ 2 4], is the ratio of the observed flux to the 
expected photospheric flux at 24 fxm. 

Using the stars identified as members with infrared excess in Figure El we plotted 
in Figure M the fraction of stars bearing disks versus spectral type for 25 Orionis (dotted 
line) and OBlb (dashed line). Error bars represent the statistical VW errors in our derived 
frequencies. Previous studies have indicated that the frequency of disks is strongly dependent 



and M) than in stars with higher masses f 


Lada & Lac 


a 


19951: 


Sicilia-Aeuilar et al. 


2005; 


Hartmann 


2005; 


Hernandez et al. 


2005; 


Carpenter et al. 


2006|; 


Briceno et al. 


2007a 


). Figure 



M also suggests that the disk frequency declines toward lower masses (spectral types later 
than Ml) showing a ma ximum arou n d K7 -M1 stars in agreement with results for the 2- 
3 Myr cluster IC 348 by lLada et al.l (120061 ). However, given the degree of uncertainty in 
each individual point in Figure El caused by the small number of disk bearing stars in each 
spectral type bin, plus the observational bias introduced by the limiting magnitude in the 
MIPS photometry, this result is not conclusive and additional data is necessary to confirm 
this trend. 
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Using Figure El the frequencies calculated for stars with evolved disks related to the total 
number of stars bearing disks in the a Orionis cluster (15.6±4.0 %), in Ori OBlb (27.8±12.4 
%) and in 25 Orionis (40±20%), indicate a clear trend toward more evolved disks in older 
stellar groups. The transitional disks candidates in these stellar groups are less frequent 
than evolved disks (8±3% for the a Orionis cluster, 17±10 % for OBlb and 10±10% for 
25 Orionis), suggesting that the transitional disk phase is relatively fast or represents an 
independent and special stage in the disk clearing period of stars. 

Disk frequencies of the sample members can be used to estimate the contamination 
level of non-members in the photometric candidates samples assuming that the photometric 
candidates bearing disks are actual members of young stellar groups. Since the overall disk 
frequencies of low mass members (6.1 ± 2.1 %) and low mass photometric candidates (7.3 
± 4.2 %) in 25 Orionis are similar, we can expect that photometric candidates sample have 
very low non-member contamination. Since the disk frequency calculated in OBlb for the 
photometric candidates (6.2 ± 3.1%) is lower than the disk frequency of members ( 13.1 
± 3.5), we can expect that around 50% of the photometric candidates are not members of 
OBlb. 



3.4. Disk evolution 



The frequency of stars bearing disks in different stellar groups is a fun ction of age, 



indicating a time scale for disk dissipation in low mass stars of 5-7 Myr (e.g., lHaisch et al. 



200ll ; IBriceno et al.ll2007al : iHernandez et al.ll2005l . 120071 ). The results presented in §3.31 are 
in agreement with this trend. In addition to this decline in disk frequency, the amount of 
infrared disk emission also decreases with age. The top panel of Figure dU] shows the median 
SED slope derived from the color [3.6] — [8.0] fo r disk bearing stars in stellar groups ranging 
in age from ~1 to ~10 Myr: 1-2 Mv r fTaurus: IHartmann et al.l l2005r). 1-3 Myr (NGC 2264 



Young et al.ll2004h . 2 -3 Myr (IC348: lLada et a 



this w ork), 7-10 Myr (25 Orionis ; IBriceno et al. 



20071 ). 4 Myr (Tr 37; ISicilia-Aguilar et al J 12009). 5 M vr (OBlb; IBriceno et al.ll2007bl . and 



2006J), 3 Myr ( a Orionis; 



Hernandez et al. 



2007bl . and this work) and 10-12 Myr (NGC 



7160; ISicilia-Aguilar et al.ll2006l ). We estimated the photospheric level using the STAR- 
PET Spitzer tool for the star templates with spectral types between K5 to M5. Error bars 
represent the quartiles given the median value for each stellar group. The median SED 
slope decreases with age, indicating a reduction of disk emission the inner regions of the 
disk. In the bottom panel, we display two sets of models described in §3.21 with different 
disk orientations, i =30 deg (dashed lines) and i=60 deg (dotted lines). The accretion rate is 



a function of age (IHartmann et al.lll998l ; iMuzerolle et al.ll2000l ; ICalvet et al.ll2005al ). so we 
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assumed M=10~ 8 Mf7) for 1 Myr o ld stellar groups, which is the mode value for accreting stars 
in Taurus (IHartmann et al.lll998l ) and M=10 _9 M Q for 10 My r old stellar groups, which is the 
mode value for accreting stars in the OB la sub- association (ICalvet et al.ll2005al ). This plot 
suggests that the expected decrease in accretion rate is not sufficient to explain the decrease 
of emission in the inner part of the disk, and it is necessary to increase the degree of settling 
(e<0.001) in the inner disk to explain the observed slopes in stellar groups with ages 5 Myr 
or older. 

Since the disk frequency and the infrared disk emission decrease with age, a correlation 
between these values can be expected. Figure [IT] confirms this suggestion, showing that 
the disk frequency is correlated with the median SED slope in the IRAC bands (correlation 
coefficient p=0.76). This plot clearly demonstrates that the disappearance of inner disks is 
related to the decrease of optical depth in th e inner disk due to the increase of dust settling 



and the decrease of mass accretion rate (see iD'Alessio et al.ll2006l ). 



Conclusions 



We have used the IRAC and MIPS data from Spitzer to study the disk frequencies and 
properties of disks around confirmed members of the Orion OBI association, 115 belong- 
ing to the ~ 10 Myr old 25 Orionis aggregate, and 106 in a region within the ~ 5 Myr 
old Ori OBlb sub association, near the Orion belt star e Ori. Using optical-2MASS color 
magnitude diagrams 41 stars were selected as additional photometric candidates of the 25 
Orionis aggregate, and 65 as additional photometric candidates of the Ori OBlb field. We 
use IRAC-MIPS diagrams to detect disk bearing stars, and to classify them as either having 
no detectable disk emission (class III), as systems with optically-thick disks (class II), or as 
objects in an intermediate phase between class II and class III systems. These intermediate 
type objects were further grouped in two categories: "transitional disks candidates", which 
have an inner, optically t hin disk region combined with an outer, optically thick disk (e.g.; 
Calvet et al.ll2002l . l2005bl ); and "evolved disk objects", i n which there is an overall decrease 



of em ission both the inner and outer regions of the disk (ILada et al.ll2006l ; IHernandez et al. 



20071 ). It is not clear if "evolved disks" are a subsequent event after transitional disks phase, 
or represent an independent evolutionary stage. 

We found that the disk frequency in the 25 Orionis aggregate (6 %) and in the Ori OBlb 
field (13 %) is mass dependent, showing a maximum value for stars with spectral type M0, 
and suggesting a decrease in the disk frequency toward higher and lo wer masses. The trend 



toward higher masses has b een observed in severa l stellar groups (e.g. IHernandez et al.ll2005 



Sicilia-Aguilar et al.l 120061 ; ICarpenter et al.l 120061 ) . The decrease in disk frequency toward 
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lowe r masses is not c onclusive, but is consistent with the results for other regions like IC 



348 (lLada et al.ll2006l ). We find that objects with evolved disks are more frequent in older 
stellar groups, while the transitional disk candidates represent a relative small fraction of the 
disk bearing stars in the various stellar groups, suggesting that the transitional disk phase 
is relatively fast. 

Comparing the disk emission in the IRAC and MIPS bands for Taurus, the a Orionis 
cluster, the Ori OB lb field and the 25 Orionis aggregate, we find that disk e mission decreases 



faster in the innermost regions of the disk; comparison with disk models from lD'Alessio et al. 



(120061 ) support this scenario. Finally, comparing the disk emission in the IRAC spectral range 
of several stellar groups ranging in age from ~1 to ~12 Myr , we find that inner disk emission 
decreases systematically with age, showing a correlation between disk frequencies and inner 
disk emission. Comparison with models using a typical accretion rate for 1 and 10 Myr 
suggests that viscous evolution alone is not sufficient to explain the decrease in the inner 
disk emission, and that large degrees of dust settling (e < 0.001) are necessary to explain 
the observed SEDs at ages 5 Myr or older. 
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Table 1. Members of the 25 Orionis aggregate 



OBla 


2MASS 


RA(2000) 


DBC(2000) 


[3.6] 


[4.5] 


[5.8] 


[8.0] 


[24.0] 


Ref 


Disk 


ID 


ID 


dcg 


dcg 


mag 


mag 


mag 


mag 


mag 


opt 


types 


9 


05224654+0134010 


80.69393 


1.56697 


12.419 ± 0.030 


12.629 ± 0.131 


12.275 ± 0.036 


12.458 ± 0.042 


. . . ± ... 


3 


cm 


25 


05224842+0140438 


80.70176 


1.67885 


12.817 ± 0.030 


12.777 ± 0.031 


12.753 ± 0.041 


12.824 ± 0.057 


. . . ± ... 


3 


cm 


47 


05225186+0145132 


80.71609 


1.75367 


13.269 ± 0.031 


13.262 ± 0.032 


13.139 ± 0.050 


13.450 ± 0.112 


... ± ... 


3 


cm 


53 


05225304+0152151 


80.72102 


1.87088 


12.416 ± 0.030 


12.348 ± 0.031 


12.240 ± 0.036 


12.258 ± 0.043 


... ± ... 


3 


cm 


905 


05245885+0125183 


81.24523 


1.42177 


12.904 ± 0.031 


12.724 ± 0.031 


12.456 ± 0.037 


11.889 ± 0.040 


9.20 ± 0.04 


3 


EV 


930 


05250192+0134563 


81.25801 


1.58232 


11.720 ± 0.030 


11.865 ± 0.030 


11.562 ± 0.032 


11.825 ± 0.036 


... ± ... 


2 


cm 


931 


05250205+0137210 


81.25855 


1.62252 


11.435 ± 0.030 


11.228 ± 0.030 


11.326 ± 0.032 


11.137 ± 0.033 


... ± ... 


1 


disk[8]? 


948 


05250362+0144121 


81.26511 


1.73670 


11.949 ± 0.030 


12.061 ± 0.030 


11.831 ± 0.034 


11.876 ± 0.036 


... ± ... 


2 


cm 



Note. — Tablc^is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown here for guidance regarding its form and content. 
References in column 10: l lBriceno et al.| j2005l ); 2 |Briceno et al.| j2007b|) ; 3 iBricefio et al.| i2007d ) 

Disk types: CIII:disk less stars; CILoptically thick disks; EV:cvolvcd disks: TD: transitional disk candidates: disk[8] ? :cxccss at 8 /^mbut no MIPS detections. 



Table 2. Members of the OBlb field 



OBlb 


2MASS 


RA(2000) 


DEC(2000) 


[3.6] 


[4.5] 


[5.8] 


[8.0] 


[24.0] 


Rcf 


Disk 


ID 


ID 


dcg 


dcg 


mag 


mag 


mag 


mag 


mag 


opt 


type 


21 


05290540-0127500 


82.27250 


-1.46390 


11.830 ± 0.030 


11.865 ± 0.030 


11.787 ± 0.034 


11.846 ± 0.038 


... ± ... 


3 


cm 


31 


05290635-0152122 


82.27648 


-1.87008 


12.751 ± 0.030 


12.707 ± 0.031 


12.632 ± 0.040 


12.595 ± 0.057 


... ± ... 


3 


cm 


63 


05290882-0125393 


82.28679 


-1.42760 


11.268 ± 0.030 


11.291 ± 0.030 


11.221 ± 0.032 


11.076 ± 0.033 


... ± ... 


1 


cm 


70 


05290925-0121227 


82.28856 


-1.35633 


12.383 ± 0.030 


12.291 ± 0.031 


12.122 ± 0.035 


11.813 ± 0.038 


... ± ... 


3 


cm 


78 


05291078-0117281 


82.29495 


-1.29115 


12.004 ± 0.030 


11.989 ± 0.030 


11.906 ± 0.034 


11.974 ± 0.040 


... ± ... 


3 


cm 


89 


05291202-0112236 


82.30010 


-1.20657 


13.138 ± 0.031 


13.085 ± 0.031 


13.120 ± 0.052 


13.067 ± 0.089 


... ± ... 


3 


cm 


14S 


05291821-0204066 


82.32590 


-2.06852 


11.267 ± 0.030 


11.152 ± 0.030 


11.058 ± 0.032 


10.917 ± 0.034 


... ± ... 


3 


cm 


209 


05292326-0125153 


82.34693 


-1.42092 


9.457 ± 0.030 


9.238 ± 0.030 


8.809 ± 0.030 


8.392 ± 0.030 


5.23 ± 0.03 


1 


CII 



Note. — Tablc[2]is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown here for guidance regarding its form and content. 
References in column 10: l lBriceno et al.| j2005l ); 2 |Briceno et al.| j2007bfi ; 3 iBricefio et al.| i2007d ) 

Disk types: CIII:disk less stars; CII:optically thick disks; EV:cvolvcd disks; TD: transitional disk candidates: disk[8] ? :cxccss at 8 /^mbut no MIPS detections. 



Table 3. Photometric candidates of the 25 Orionis aggregate 



OBla 


2MASS 


RA(2000) 


DEC(2000) 


[3.6] 


[4.5] 


[5 


.8] 


[8.0] 


[24.0] 


V 


\'I 


Disk 


ID 


ID 


dcg 


dcg 


mag 


mag 


mag 


mag 


mag 


mag 


mag 


type 


146 


05230905+0125355 


80.78774 


1.42654 


11.383 ± 0.030 


11.355 ± 0.030 


11.278 


± 0.032 


11.284 ± 0.034 


. . . ± ... 


15.75 ± 0.04 


2.43 ± 0.05 


cm 


297 


05233109+0144079 


80.87958 


1.73555 


11.005 ± 0.030 


11.079 ± 0.030 


11.028 


± 0.031 


10.967 ± 0.032 


. . . ± ... 


13.91 ± 0.03 


1.06 ± 0.05 


cm 


359 


05234182+0152261 


80.92428 


1.87394 


13.144 ± 0.031 


13.088 ± 0.031 


12.986 


± 0.041 


13.039 ± 0.055 


... ± ... 


19.11 ± 0.05 


3.09 ± 0.09 


cm 


427 


05235215+0136314 


80.96732 


1.60873 


11.103 ± 0.030 


11.191 ± 0.030 


11.155 


± 0.032 


11.074 ± 0.032 


... ± ... 


13.84 ± 0.03 


0.98 ± 0.05 


cm 


458 


05235854+0151255 


80.99394 


1.85709 


12.114 ± 0.030 


12.047 ± 0.030 


12.003 


± 0.034 


12.043 ± 0.040 


. . . ± ... 


17.29 ± 0.04 


2.58 ± 0.07 


cm 


477 


05240118+0128236 


81.00493 


1.47324 


12.316 ± 0.030 


12.289 ± 0.031 


12.123 


± 0.035 


12.165 ± 0.041 


. . . ± ... 


16.87 ± 0.04 


2.49 ± 0.06 


cm 


543 


05241034+0155024 


81.04312 


1.91735 


10.322 ± 0.030 


10.428 ± 0.030 


10.365 


± 0.031 


10.309 ± 0.031 


. . . ± ... 


13.16 ± 0.03 


0.80 ± 0.05 


cm 


1626 


05265473+0144337 


81.72806 


1.74271 


10.168 ± 0.030 


10.273 ± 0.030 


10.061 


± 0.031 


10.069 ± 0.031 


9.83 ± 0.06 


13.43 ± 0.03 


0.86 ± 0.05 


EV 



Note. — Tablc[3]is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown here for guidance regarding its form and content. 
Disk types: CIILdisk less stars; CIPoptically thick disks: EV:cvolvcd disks; TD: transitional disk candidates; disk[8] ?:cxccss at 8 /xmbut no MIPS detections. 



Table 4. Photometric candidates of the OBlb field 



OBlb 2MASS RA(2000) DEC(2000) [3.6] [4.5] [5.8] [8.0] [24.0] V VI Disk 

ID ID deg deg mag mag mag mag mag mag mag type 



74 


05290984- 


0208250 


82.29101 


-2.14030 


12.269 


± 


0.030 


12.217 


± 


0.031 


12.141 


± 


0.036 


12.113 


± 


0.042 


. .. ± . . 


18.26 ± 0.05 


2.89 ± 0.09 


cm 


187 


05292160- 


0201546 


82.34002 


-2.03184 


11.142 


± 


0.030 


11.121 


± 


0.030 


11.073 


± 


0.031 


11.025 


± 


0.034 


. . . ± . . 


15.91 ± 0.06 


2.04 ± 0.09 


cm 


203 


05292300- 


0126559 


82.34584 


-1.44888 


13.500 


± 


0.131 


13.262 


± 


0.032 


13.186 


± 


0.042 


13.184 


± 


0.085 


. .. ± . . 


18.89 ± 0.05 


2.83 ± 0.08 


cm 


207 


05292313- 


0149203 


82.34638 


-1.82231 


13.109 


± 


0.031 


13.105 


± 


0.031 


13.059 


± 


0.051 


13.451 


± 


0.153 


. . . ± . . 


18.83 ± 0.05 


2.63 ± 0.08 


cm 


222 


05292426- 


0207354 


82.35112 


-2.12651 


12.953 


± 


0.031 


13.055 


± 


0.031 


12.730 


± 


0.042 


12.825 


± 


0.063 


. . . ± . . 


17.78 ± 0.05 


2.27 ± 0.08 


cm 


236 


05292591- 


0144580 


82.35797 


-1.74945 


10.810 


± 


0.030 


10.832 


± 


0.030 


10.757 


± 


0.031 


10.707 


± 


0.032 


. . . ± . . 


14.64 ± 0.04 


1.41 ± 0.06 


cm 


276 


05293010- 


0114446 


82.37544 


-1.24573 


11.116 


± 


0.030 


11.217 


± 


0.030 


11.175 


± 


0.032 


11.107 


± 


0.034 


. . . ± . . 


14.37 ± 0.00 


1.25 ± 0.00 


cm 


283 


05293049- 


0121500 


82.37707 


-1.36389 


11.303 


± 


0.030 


11.376 


± 


0.030 


11.314 


± 


0.032 


11.309 


± 


0.036 


. . . ± 


14.56 ± 0.04 


1.30 ± 0.06 


cm 


342 


05293578- 


0148046 


82.39909 


-1.80129 


11.138 


± 


0.030 


11.185 


± 


0.030 


11.129 


± 


0.031 


11.144 


± 


0.034 











Note. — Tablc^is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown here for guidance regarding its form and content. 
Disk types: CIILdisk less stars; CII:optically thick disks; EV:evolvcd disks; TD: transitional disk candidates; disk[8]?:exccss at 8 pmbut no MIPS detections. 
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Fig. 1. — IRAC image of the 25 Orionis aggregate. This is a 3-color composite of IRAC 
images (3.6 (blue), 4.5 (green), and 8.0 (red) fim) illustrating the space distribution of 
members in this ag gregate. Circles represent stars confirmed as low mass members using 
spectroscopic data ( Bricefio et al. 2005 . 2007bl lc1) . squares represent photometric low mass 
candidates of 25 Orionis. Solid symbols indicate stars bearing disks (see §3.ip . Intermediate 
mass members selected in Hernandez et al. J2006h are represented as diamonds, open symbols 
represent stars without disks and soli d symbols represent deb ris disk candidates. The star 
indicates the Herbig Ae star V346 Ori (jrlernandez et al" 2006 ) and the triangle indicates the 
Be star 25 Ori. [See electronic edition of ApJ for a color version of this figure.] 
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Fig. 2. — IRAC image of the OBlb region. This is a 3-color composite of IRAC images, 3.6 
(blue), 4.5 (green), and 8.0 (red) (im. Symbols a re as in Figured! The star indicates the 
Herbig Ae star HD290543 Jliernandez et alJl2006h [S ee electronic edition of ApJ for a color 
version of this figure.] 



-24- 




- > 



2 3 4 5 6 
[V-J] 



14 
16 
18 

20 



I I |\ I M\l | I I I I | I I I I | I I I I | I I I I 

\V* OB lb J 




\ \ 

I I I I I I I I I I I I I I I I I Nl I I I I I I l N 



14 
16 
18 

20 



25 Ori _I 




14 
16 
18 

20 



i \ i \ k i — i — i — i — I — i — i — i — i — I — i — i — r 

h \ \ \ 

\\ OBlb J 




2 3 4 5 6 
[V-J] 



Fig. 3. — Color magnitude diagrams illustrating the selection of low mass photometric can- 
didates in 25 Orionis (upper panels) and in OBlb (lower panels). Photometric candidates 
(open squares) have similar distribution than m embers (open circles) confirmed using spec- 
troscopic methods ( Bricefio et al. 2005 . 2007bl jcl). By co mparison, the ZA MS (solid lines) 
and the isochrones at 10 Myr and 5 Myr ( dashed lines; ISiess et al.l I200CH ) ar e displayed at 
the d i stance of each stellar group; 3 30 pc for 25 Orionis and 440 pc for OBlb (jBriceno et al. 
20051 . l2007bl ; iHernandez et al.l 120051 ) . Long dashed lines represent the median colors for the 
spectroscopic members. [See electronic edition of ApJ for a color version of this figure.] 
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Fig. 4. — Diagrams illustrating the detection of sources with infrared excess above the 
photospheric levels (solid lines) in 25 Orionis. The top panel shows the IRAC SED slope 
diagram indicating stars with excess at 8 ^m. The bottom-left panel shows a 2MASS-MIPS 
color-color diagram indicating stars with excess emission at the MIPS 24 /im band. The 
bottom-right panel shows an IRAC color-color diagram indicating objects with infrared ex- 
cess at [4.5] and [8.0]. Members and photometric candidates are displayed as circles and 
squares, respectively. Large circles and large squares represent stars with infrared excess 
at 24 /im. The solid histograms repr esent the distribution of IRAC SED slopes for stars 
with optically thick di sks in Taurus (IHartmann et al.l 120051 ) and in the a Orionis cluster 
([Hernandez et al.l 120071 ); dashed lines represent the class II region limit based on these his- 
tograms. The open histograms represent the distribution of stars with no IRAC excesses. 
The dashed histogram is the distribution of stars with evolved disks in the a Orionis cluster 
([Hernandez et al.l 120071 ) [See electronic edition of ApJ for a color version of this figure.] 




Fig. 5. — Diagrams illustrating the detection of sources with infrared excess above the 
photospheric levels (solid lines) in the OBlb field. Symbols are similar to Figure HI [See 
electronic edition of ApJ for a color version of this figure.] 
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Fig. 6. — SED slope K— [5.8] versus SED slope K— [24] for members (circles) and photometric 
candidates (squares) bearing disks in 25 Orionis (left panel) and in OBlb (right panel). Error 
bars represent the quartiles of the stellar groups, Taurus (1 Myr), the a Orionis cluster (3 
Myr), OBlb (5 Myr), and 25 Orionis (10 Myr). Stars repre s ent the transitional disk objects, 



Coku Tau4 TW Hya, and GM Aur (ID'Alessio et al.ll2005al : llJchida et al.ll2004l : ICalvet et al. 



20021 . |2005bJ) . Photospheric regions derived for K5-M5 stars using STAR- PET Spitzer tool 



are re presented as dashed areas. Disk bearing stars in the a Orionis cluster ([Hernandez et al. 



20071 ) are represented with the symbol X. The horizontal dotted line represents the lower 
quartile of the stars bearing disks in the a Orionis cluster, ~96 % of stars with disks in 
Taurus are located above this line; we define the area above this line as the class II region 
(CII in Table 1-4). The vertical dotted line represents the lower quartile of the stars bearing 
disks in the a Orionis cluster. We use this limit to separate transitional disk candidate (TD 
in Table 1-4) from stars with evolved disks (EV in Table 1-4). 
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Fig. 7. — Spectral energy distributions for selected disk bearing stars, illustrating the diver- 
sity of disks fou nd at 5 and 10 Myr. Do tted line represents the median SEDs of disk bearing 
stars in Taurus ( IHartmann et al.ll2005l ). Dashed line represents the median photosphere of 
stars in the spectral type range K5-M5 (IKenyon &: Hartmannlll995l ). We display examples 
of the disk types found in Figure [6j stars with optically thick disks (CII), transitional disk 
candidates (TD) and evolved disks objects (EV).[See electronic edition of ApJ for a color 
version of this figure.] 
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Fig. 8. — Theoretical SED slopes as function of the degree of settling. Each panel shows 
three curves corresponding to different accretion rates (M=10 _7 ,M=10 -8 and M=10 -9 ). The 
parameter e represents the degree of settling in each model. The central object is a typical 
T Ta uri star with M=O.6M , Age=l Myr, corresponding to a spectral type K7 (ISiess et al. 
20001 ) . Left and middle panels represent disks oriented with an angle of 30deg and 60 
deg, respectively. In the right panels, we display the qua rtiles of the disk population in 
Taurus (1-2 Myr; lHartmann et al.ll2005l ; iFurlan et al.ll2006l ). the a Orionis cluster (~3 Myr 
Hernandez et al.ll2007l ). and the stellar groups studied in this contribution, OBlb (~ 5 Myr) 
and 25 Orionis (~ 10 Myr). [See electronic edition of ApJ for a color version of this figure.] 
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Fig. 9. — Disk frequency as function of spectral type for the 25 Orionis aggregate (dotted 
line) and the OBlb field (dashed line). The error bars represent the Poissonian statistical 
uncertainties 
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Fig. 10.— Medi an SEP slope in t 



10 



Age (Myr) 



le IRA C spectral ra nge versus stellar a ges of diff erent stellar 



group s: Taurus (Hartm ann et al 
200fij l. a Orionis (jflernandez et al. 



and 25 Orionis and OBlb ( Briceno et al. 



2005), NGC 2264 flYoung et al.ll2004h. IC 348 flLada et al. 



20070. Tr37 and NGC 7160 (ISicilia-Aguilar et alJl2006h 



2007bl . and this work) . The dashed region indicates 



the photospheric levels estimated using STAR-PET Spitzer tool. It is apparent that the 
infrared disk emission decreases with the age of the stellar groups. The lower panel shows 
a comparison with theoretical values described in §3.21 Two set of models are displayed 



showing different orientation of the disks, 30 deg (dashed lines) and 60 deg (dotted lines). 
We assume canonical values for accretion rate at 1 Myr and 10 Myr (see text). Different 
lines describe different degree of settling (e=l, 0.1, 0.01 and 0.001) showing flatter slopes for 
models with less degree of settling. [See electronic edition of ApJ for a color version of this 
figure.] 
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Fig. 11. — IRAC disk Frequency versus the median SED slope [3.6] — [8.0] for different stellar 
groups (see Figure fTDl . These measurements are correlated (correlation coefficient p=0.76) 
indicating more flaring disks in stellar groups with higher disk frequency. 



